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We present a statistical analysis on the performance of a protocol for the faithful transfer of a 
quantum state in finite qubit or spin chains, in the presence of diagonal and off-diagonal disorder. It 
is shown that the average-state fidelity, typically employed in the literature for the quantification of 
the transfer, may overestimate considerably the performance of the protocol in a single realization, 
leading to faulty conclusions about the success of the transfer. 
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I. INTRODUCTION 

The faithful transfer of quantum states between two 
or more spatially separated qubits of a quantum network 
has attracted considerable interest over the last decade 
In its simplest form, the problem of quantum state 
transfer (QST) pertains to the quest for Hamiltonians 
(protocols) that ensure the transfer of a state between 
the two ends of a qubit chain at a prescribed time Q. 
The qubit chain may be represented by an actual spin 
chain in a liquid or solid-state NMR system [E |j| , ni- 
trogen vacancies in diamond [5| , or by other realizations 
pertaining to optical lattices |6|, coupled quantum dots 
@, 111 i superconducting qubits [E EE| , or waveguides [ll| ■ 
The first proof-of-principlc experiment on a QST Hamil- 
tonian with engineered couplings has been presented re- 
cently in the framework of waveguides (llT | . 

Disorder is expected to be present in any physical re- 
alization of qubit (spin) chains, irrespective of the ex- 
perimental platform. Thus, theoretical investigations on 
the robustness of QST Hamiltonians against disorder is 
an essential step toward any experimental realization, as 
well as for the classification of different protocols, with re- 
spect to their performance under realistic conditions. To 
the best of our knowledge, so far such investigations have 
relied on an input-independent fidelity obtained by aver- 
aging over all possible input states [EEE Il2| - [l5| . This 
average-state quantity is evaluated for a large number 
of independent realizations (pertaining to randomly cho- 
sen disorder for the diagonal and off-diagonal elements 
of the Hamiltonian), and all the analysis of the perfor- 
mance of a protocol is based on the ensemble averaged 
input-independent fidelity. 

In order for a QST Hamiltonian to be reliable and use- 
ful, however, its performance for a certain level of disor- 
der has to be within the acceptable levels for every single 
realization, irrespective of the input state. The disorder 
as well as the qubit state to be transferred may vary from 
realization to realization, but the efficiency and the suc- 
cess of the protocol under consideration have to be guar- 
anteed. Our aim in this work is to investigate whether 
and to what extent the average-state input-independent 
fidelity that has been used extensively in the literature, is 
capable of describing the performance of a Hamiltonian 



in a single realization. To this end, we consider one of 
the most studied QST Hamiltonians 0, @, El EE BSl , 
the robustness of which in the presence of disorder has 
been also investigated in various contexts 0, EE EH~EI| • 
In Sec. |H] we outline our model, whereas Sees. IIIII 
and IIVI are devoted to a thorough statistical analysis on 
the performance of the protocol in a single realization. 
It is shown that for a large class of input states the en- 
semble averaged input-independent fidelity overestimates 
the performance of the protocol, and may lead to faulty 
conclusions with respect to the success (or failure) of the 
transfer. Alternative forms of the fidelity turn out to 
be more reliable measures for the quantification of QST. 
Our conclusions are summarized in Sec. IV] 



II. MODEL 

Our analysis pertains to QST spin-chain Hamiltonians 
of the form 
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where N the length of the spin chain, o^ v ' z are the Pauli 
spin operators for the fcth spin, Ek is the energy separa- 
tion between the spin-up and spin-down states playing 
the role of the local "magnetic field" , and Jk is the time- 
independent nearest-neighbor spin-spin interaction. This 
is a spin-chain Hamiltonian of the XX type, which is iso- 
morphic to the Hubbard Hamiltonian for non-interacting 
spinless fermions or hard-core bosons plj 
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where aj. (a^) is the particle creation (annihilation) op- 
erator at fcth site with energy and Jk now plays the 
role of tunnel coupling between adjacent sites k and 
From now on, the qubit basis states for a single spin (par- 
ticle) are denoted by {|0), |1)}. 

Our task is the faithful — ideally perfect — transfer 
of an input qubit state from the first to the iVth site of a 
chain that is governed by a Hamiltonian of the form ([1]) 
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FIG. 1. (Color online) Probability distribution of the fidelity for N = 12 and 1000 independent realizations. The histograms (a- 
e) correspond to F^(\f3\ 2 ,p,Aip) for |/3| 2 = 1.0, 0.8, 0.6, 0.5, and 0.4, respectively. The histogram (f) corresponds to F(p,Aip). 
Other parameters: a v = = 0.1. 



[or ([2])]. at a prescribed time. Various Hamiltonians have 
been found to achieve this task [l[, one of which will be 
adopted later on for our purposes. In general, the input 
state need not be a pure state, and a measure for the 
fidelity of the transfer is [22|, HH 

F(p,a)= (TrvV/2^1/2) 2 , (3) 

where a is the input state, and p is the state at the out- 
put. The minimum fidelity 

F min = imn{F{p,a)}, (4) 

is obtained by taking the minimum over all the possi- 
ble input qubit states. Using the joint concavity of the 
fidelity one can show that it is sufficient to take the min- 
imum over all the possible pure qubit states [22|, [23[ i.e., 

F m in = min{F^} (5a) 

•0 

where 

F 4 , = Mp|V>, (5b) 

and 

|V) = a|0)+/3|l), (6) 

with \a\ 2 + \/3\ 2 = 1 and a, j3 £ C. 

Assume that the entire chain is initially prepared in 
the ground (vacuum) state, and the first spin (site) is 
prepared in state \tp). The Hamiltonian ([T]) [or ©] pre- 
serves the number of spin- [or particle-] excitations, and 
thus we need to consider only the zero |0) = FJ J=1 |0)j 
and single excitation \j) = &j~\0) (d]|0)) sectors of 



the total Hilbert space. Then the entire chain is ini- 
tially in state 1^(0)) = a\0) + f3\l) and evolves in 
time as |*(t)) = #(t)|¥(0)) = a|0) + /?E^i A Ml), 
where & (t) = exp [A J4?t\ is the evolution operator, and 

Apparently, only the states in the single excitation sec- 
tor evolve in time, while the ground (or vacuum) 
state |0) remains unchanged. We are interested in the 
transfer of the input state from the first to the ATth site. 
The reduced operator for the Nth. site at time t is 

p(t) = TV[|*(t)X*(t)|] (7) 
= (1-|/3| 2 |^| 2 )|0)(0| + |/3| 2 |^| 2 |1)(1| 

+a/3*A* N \Q)(l\+a*pA N \l)(0\, (8) 

where 

A N (t)={N\^(t)\l) = y/plfje^, (9) 

and pit) is the probability for the excitation transfer. 
Using Eqs. (|5b]). (JSJ), and © we obtain 

F^(|/3| 2 ,p, ip) = 1 + |/3| 2 [-1 - p + 2>cos(^)] 

+ [2 13 -2^cos(^)]|/3| 4 . (10) 

This is the fidelity in a single realization of the protocol, 
for a given input state \ip). Clearly, it depends on the 
input state (through /3), the chosen QST Hamiltonian 
and the presence of imperfections through p and (p. 

In the absence of disorder, at the end of an ideal single 
realization of a QST protocol, one would have a well- 
defined amplitude An = y/poe llpa , with p < 1 (equality 
holds for perfect-state-transfer protocols), whereas the 
phase is fixed and known, and thus, in principle, cor- 
rectable. Without loss of generality, we can assume that 
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one compensates for ipo at the end of any (ideal or non- 
ideal) realization of the transfer, and thus the fidelity 
(ITUl) can be rewritten as El-fll 



1 + 1/ 
+ [2p 



| 2 [-l-p + 2^cos(A^)] 
-2^cos(A^)]|/3| 4 , (11) 



where Aip = ip — ipo. Apparently, for an ideal realization 
of any faithful QST protocol one has p = po and Aip = 0. 

Any physical implementation of a QST protocol will in- 
evitably deviate from the ideal scenario due to the pres- 
ence of disorder. In principle, one can distinguish be- 
tween static and dynamic disorder, which manifest them- 
selves through the randomization of the diagonal (ener- 
gies) and off-diagonal (couplings) elements of the Hamil- 
tonian. In the case of dynamic disorder, such a random- 
ization varies within the time interval that the transfer 
takes place. On the other hand, in the case of static 
disorder the randomization remains practically constant 
during a single realization of the transfer. The origin 
of such a disorder can be attributed to manufacturing 
errors, or to dynamic unpredictable factors (e.g., tem- 
perature fluctuations, stray fields, etc), which introduce 
a time-dependent disorder that does not vary apprecia- 
bly over the time scale of a single transfer. It varies, 
however, appreciably from realization to realization in a 
random way, and thus one cannot compensate for the 
unpredictable changes. 

The main effect of the disorder, static or dynamic, is 
to randomize p, ip (around po and ipo > respectively) , and 
thus the single-realization fidelity F^p also becomes a ran- 
dom number that depends on the input state \ip). For the 
rest of this section we will focus on static disorder. To 
the best of our knowledge, so far analogous studies in the 
literature have relied on the input-independent fidelity 
ilii 



F(p, A V ) = i + | 



v ^cos(A<^) 



(12) 



obtained by averaging ([lip over all the possible input 
states. This average-state fidelity still varies from real- 
ization to realization and by taking an ensemble average 
over many realizations, the performance of a particular 
QST protocol is usually quantified by the ensemble aver- 
aged input-independent quantity (F(j>, Aip)}. 

A faithful QST Hamiltonian should operate reliably 
in every single realization and irrespective of the input 
state. A question therefore is how reliably F(p, Aip) de- 
scribes a single realization, where one typically has the 
transfer of a particular, and a priori unknown, qubit 
state \ip). To answer this question, one has to compare 
F(p, A(p) to the single-realization fidelity F^,(\(3\ 2 ,p, Aip) 
for a particular QST Hamiltonian. In the following sec- 
tion we present such an analysis. 
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FIG. 2. (Color online) The average fidelities (F^(|^| 2 ,p, Aip)) 
(solid curves) and (F(p, Atp)\ (horizontal solid lines) as func- 
tions of \/3\ 2 , for different values of TV: (a) N = 12; (b) N = 18; 
(c) N = 25; (d) N — 31. The averages have been obtained 
on 1000 independent realizations and the error bars show the 
standard deviations as obtained from the corresponding dis- 
tributions of F(p,Aip) and F^(\/3\ 2 ,p, Aip) for various values 
of |/3| 2 . The standard deviation of F(p, Aip) is independent of 
|/3| 2 . The dashed horizontal line shows the classical threshold 
F c i = 2/3. Other parameters as in Fig. [Tj 



III. FIDELITY STATISTICS 

The previous formalism applies to any QST Hamilto- 
nian of the form (TTJ [or @]. For the sake of concrete- 
ness, however, the following analysis is focused on a well- 
studied QST Hamiltonian involving 0, 0, M, 



e k = s, and J k = J \Jk(N - k), 



(13) 



(see also Ref. [l_l| for the first related proof-of-principle 
experiment). It has been shown that this Hamiltonian 
under ideal conditions ensures the perfect transfer of a 
qubit state between the two ends of the chain, at time 
t = 7r/(2 J ) i-e., 



A n (t) 



H) 



N-l 



so that po = 1 [241 ]. In any physical realization of such a 
Hamiltonian, it is expected to exist an upper bound on 
the achievable coupling strengths J max = max{Jfe}, im- 
posed by technological or physical constraints. Through- 
out our simulations we worked with the dimensionless 
quantities e k = fffc/Jmax and J k = J k /J ma , x . We per- 
formed a large number of realizations, each one pertain- 
ing to a sequence of independent random variables r\ k and 
£fc, normally distributed around zero, and with standard 
deviations a v and a^ 7 respectively. The energies and the 
couplings entering our QST Hamiltonian were random- 
ized as follows 



e k ^e + ij kl and J k -> J k {l + £ fc ). 



(14) 



In each realization, we kept track of An (t) obtaining thus 
a large sample of the parameters entering the fidelities. 
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FIG. 3. (Color online) (a) The average fidelities 
(F^(\j3\ 2 ,p, Aip)) for various values of |/3| 2 (solid curves with 
symbols) and \F(p, Atp) J (dashed curve with circles) are plot- 
ted as functions of TV. (b) As in (a) for the standard devi- 
ations, (c) The intervals Aj for which the average fidelity 
(F^(\/3\ 2 ,p, Atp)} is below (F(p, Atp)) (dot-dashed curve with 
triangles) and the classical limit (solid curve with squares), 
as functions of TV. (d) The probability of failure for various 
values of |/3| 2 , as a function of TV. The vertical arrows mark 
the values of TV for which (F(p, Atp)) — F c i. The inset of 
Fig. [3jb) indicates the values of |/3| 2 that correspond to the 
depicted curves. Other parameters as in Fig. [2] 
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FIG. 4. (Color online) The average fidelities (F(p,Aip)) (a) 
and {F^\ 2 ,p,Atp)) for \/3\ 2 = 0.8 (b), \(3\ 2 = 0.6 (c) and 
\/3\ 2 = 0.4 (d), as functions of disorder for TV = 15. The aver- 
ages have been obtained on 1000 independent realizations. 



As mentioned above, the performance of the protocol 
in a single realization can be studied cither in terms of 
the average-state fidelity F(p,Ap), or the actual single- 
realization fidelity F^(|/3| 2 ,p, Atp). Due to the presence 
of disorder, both of them fluctuate from realization to 
realization, and in Fig. [T] we show the corresponding 
distributions after 1000 realizations, for various input 
states (i.e., various values of |/3| 2 ). Clearly, in many 
cases the distributions are rather broad, which means 
that the performance of the protocol in a single realiza- 
tion may deviate considerably from the mean fidelities 
(F) and (Fy). Moreover, for |/3| 2 > 0.5 the distributions 
of i^0(|/3| 2 ,p, Atp) [see Fig. [lja-c)] are wider than the 
distribution of F(p, Atp) [see Fig. [TJf )]] , and at the same 
time they are peaked at lower fidelities; the opposite is 
true for |/3| 2 < 0.5. This is a typical behaviour of the 
distributions for fixed TV, which suggests that F tends 
to overestimate (underestimate) the performance of the 
protocol in a single realization for states with \/3\ 2 > 0.5 
(1/3 1 2 /-5 0.5), respectively. 

In Fig. [2] we plot the ensemble averaged fidelities 
(F,/,(|/3| 2 ,p, Atp)) and (F(p,Atp)) as functions of |/?| 2 , to- 
gether with the classical limit F c i = 2/3 (dashed line) 
(25| , for various values of TV. In addition we plot the stan- 
dard deviations of the distributions of F^(|/3| 2 ,p, Atp) 
and F(p,Atp), for various |,$| 2 (vertical bars). In agree- 
ment with our observations on Fig. Q] we see that for 
1/3 1 2 < 0.5, (F$) is above (F) and decreases with increas- 
ing |/3| 2 , crossing (F) at |/3| 2 sa 0.5. In addition, the 



standard deviation of F^ increases with increasing |/3| 2 , 
which implies a larger dispersion of the single realiza- 
tion fidelities around the mean value (ify). As a result, 
(F(p, Atp)') is not a reliable measure of the performance 
of the protocol under consideration, since for a broad 
range of input states with \/3\ 2 > 0.5, most likely the fi- 
delity of the transfer in a single realization will be well 
below (F(p, Atp)'). More importantly, for input states 
with |/3 1 2 close to 1 and for relatively large values of TV 
(depending on the disorder under consideration) , it is also 
highly probable that the fidelity of the transfer in a single 
realization of the given QST protocol will be below the 
classical limit F c i , which implies that the particular state 
transfer has failed, yet (f 1 ) is well above F c i indicating 
faithful QST [e.g., see Figs. E^c) and (d)]. 

To make this point clearer, in Fig. |3ja) we plot the 
dependence of (F^(\/3\ 2 7 p, Atp)) and (F(p,Ap)') on the 
number of sites in the chain, for a fixed disorder. In all 
cases, the ensemble averaged fidelities decrease with TV, 
but for input states with \(3\ 2 > 0.5, (F^) decreases con- 
siderably faster than (F'), crossing the classical limit at 
smaller values of TV. Moreover, the corresponding stan- 
dard deviations increase with TV [see Fig. [3] (b)], and in 
general the ones of F^ for |/3| 2 > 0.5 are larger than the 
standard deviation of F, for all the values of TV. 

As depicted in Fig. [21 for any fixed TV, 
(F^(|/3| 2 ,p, Atp)) = 1 at |^| 2 = 0, and decreases with 
increasing |/3| 2 , whereas (F(p,A<p)} and F c i are inde- 
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FIG. 5. (Color online) The probability of failure as estimated 
by (F(p,A<p)) (a), and (F 4 ,{\P\\p, Atp)) for |/3| 2 = 0.8 (b), 
|/3| 2 = 0.6 (c) and |/3| 2 = 0.4 (d), as a function of disorder for 
TV = 15. The averages have been obtained on 1000 indepen- 
dent realizations. Successive contours differ by 0.1 and the 
color scale is as in Fig. [4] 
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FIG. 6. (Color online) The quantities Aq (a) and Ai (b) as 
functions of disorder for TV = 15. The averages have been 
obtained on 1000 independent realizations. 



pendent of |/3| 2 . The values of \/3\ 2 e [0,1] at which 
(F^(\f3\ 2 ,p, Atp)} crosses (F(p, Atp)) and F c i, are of par- 
ticular interest since they essentially characterize the 
class of input qubit states for which (F^(|/3| 2 ,p, Atp)) be- 
comes smaller than the other two fidelities. For such a 
characterization, we can introduce here the quantity 



1 - B 



(c) 



if < B^ < 1 



0, otherwise 



(15) 



where Bj, with j £ {0,1}, are the roots of the poly- 
nomials S (\(3\ 2 ) = (i^(|/3| 2 ,p, Ay.)) - (F(p,Atp)) and 
Si(\f3\ 2 ) = (F^(\j3\ 2 ,p, Atp)) - F cl , with respect to |/3| 2 . 



For all the tested parameters and noises in our simu- 
lations, these polynomials had at most one root in the 
interval [0,1], and thus Aj were uniquely defined. As 
shown in Fig. EHc), Ao does not vary considerably with 
N, and throughout our simulations it was around 0.5. 
This means that for all N, there is always a class of input 
qubit states with \/3\ 2 > B^ c \ where B$ m 0.5, for which 
{F(p, Atp)) overestimates the performance of the proto- 
col relative to (F^(\f3\ 2 ,p, Atp)). Moreover, for relatively 
small values of N, there is no root for the polynomial 
Si in [0,1], and thus Ai = 0, which means essentially 



that (F^(|/3| 2 ,p,A(^)) > F ch for all 



G [0, 1]. As we 



increase N, however, the polynomial Si acquires a root 
b[ c ' € [0, 1], which appears close to 1 for N w 23, and 
moves rapidly toward 0.4, with increasing N. As a result, 
the range of values of (3 for which (F^,(\f3\ 2 7 p, Atp)) < F c [ 
is getting wider with increasing N > 23, whereas at the 
same time, one has {F(jp, Atp)) < F c i only for N > 73. 
Hence, for 23 < N < 73 the fidelity (F(p,Atp)) pre- 
dicts incorrectly the transfer of input qubit states with 
\P\ 2 > b[ c \ wher eas the transfer of such qubit states 
under the particular QST protocol is most likely impos- 
sible (i.e., the fidelity in a single realization is below the 
classical limit). 

Of course, as mentioned before, there are statistical 
deviations from realization to realization and thus in 
our simulations we have also kept track of the num- 
ber of realizations for which F,p(\/3\ 2 ,p, Atp) > F c i and 
F(p,Atp) > F c i. Dividing by the total number of real- 
izations, we thus have an estimate of the probabilities of 
failure for the protocol under consideration and for vari- 
ous input states. As shown in Fig. E^d), for all the states 

with |/3 1 2 > B^ (see the curves with squares and trian- 
gles corresponding to \(3\ 2 = 0.8 and 0.6, respectively), 
the probability of failure that is estimated according to 
Fjf,(\/3\ 2 ,p, Atp) can be considerably higher than the one 
that is obtained from F(p, Atp) (dashed curve with cir- 
cles). The fidelity F(p,Atp) predicts reliably the prob- 
ability of failure only in two cases. Namely, either for 
states with \/3\ 2 < B^ c ' for all TV (sec curve with stars 
corresponding to |/3| 2 = 0.4), or for all the input qubit 
states when N < 10. 

Our discussion so far has been restricted to a partic- 
ular combination of diagonal and off-diagonal static dis- 
order (i.e., fixed a v and a^). In Fig. H{a), we present 
a contour plot of {F(p,Atp)) as a function of <r v and 
0£, whereas Figs. |3Jb-d) are the corresponding plots 
of (F^(|/3| 2 ,p, Atp)), for various values of |/3| 2 . In all 
cases, the fidelities drop with increasing disorder, but for 
|/3| 2 = 0.6 and |/?| 2 = 0.8 <F^(|/3| 2 ,p, Atp)) decreases 
considerably faster than (F(p,Atp)), whereas the oppo- 
site is true for |/3| 2 = 0.4. Analogous contour plots for 
the probability of failure are depicted in Fig. [5] Again 
we see that as we increase the disorder for |/3| 2 = 0.6 and 
\f3\ 2 = 0.8, the probability of failure for (F,/,(|/3| 2 ,p, Atp)) 
increases considerably faster than the probability of fail- 
ure for (F(p, Atp)). The latter seems to overestimate the 
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FIG. 7. (Color online) The values of |/3| 2 that minimize F^(\/3\ 2 ,p, Atp) (a) as well as the fidelity differences -F m in(p, Aip) — 
F(\,p, Atp) (b) and Fmi n (p, Aip) — F(p,Atp) (c), for various combinations of p and Atp. The plots are mirror-symmetric with 
respect to Atp = 0, and can be extended to Atp < by the transformation Atp — y —Atp. 



probability of failure only when |/3| 2 < 0.5 for the spe- 
cific value of N. Finally, Fig. [5] shows that the interval 
of values for \f3\ 2 where (F^(\f3\ 2 ,p, Atp)) < (F(p,Atp)) 
does not vary appreciably with the disorder (it is around 
0.5). On the contrary, the corresponding interval where 
(F^(\f3\ 2 ,p, Aip)) < F c i becomes wider with increasing 
disorder. 

To summarize, all of the above observations confirm 
the fact that the average-state fidelity F(p,ip), usually 
employed in studies on QST, is not a reliable measure 
for the performance of QST protocols in a single real- 
ization for all the possible levels of disorder, and for all 
possible input qubit states. It tends to overestimate the 
performance of QST protocols for a broad class of qubit 
states, and it may even lead to faulty conclusions with 
respect to the classical limit. As will be discussed in the 
following section, a more reliable measure for the perfor- 
mance of QST protocols is provided by F m { n7 or even by 
the probability p in certain cases. 



IV. RELIABLE MEASURES OF STATE 
TRANSFER 

The most reliable measure for the quantification of 
QST is the minimum fidelity (|5a[) . The minimization 
of the single-realization fidelity F^ , over all the possible 
input qubit states is straightforward in the absence of dis- 
order where p = po and Atp = 0. Indeed, one can readily 
show that Fjf,(\/3\ 2 ,po, 0) is a monotonically decreasing 
function of |/3| 2 G [0,1], which attains its minimum at 
|/3 1 2 = 1 and is equal to po. This means that in the 
case of ideal realizations, the transfer of the excitation 
(probability) implies transfer of the quantum state with 
precisely the same fidelity. In the case of perfect QST 
Hamiltonians, po — 1 for all input states. 

In general, the transfer of the excitation is quantified 
by the probability p, which in the presence of disorder 
static or dynamic, is a random variable in the interval 



[0,po]j whereas Aip is a random variable in [— it, it]. In 
this case, the transfer of the excitation does not nec- 
essarily imply transfer of the state (which in addition 
to bit information also carries phase information) [27j . 
For relatively weak disorder, the distributions of p and 
Aip are expected to be peaked at the corresponding ideal 
values po and 0, respectively, and how fast the distribu- 
tions drop depends on the type of noise under considera- 
tion. From the mathematical point of view, irrespective 
of the Hamiltonian and the disorder under consideration, 
F,f,(\f3\ 2 ,p, Atp) is a stochastic function, which for any 
given p <G [0, 1] and Aip <G [— 7r, tt] attains its maximum 
for (3 = 0. Moreover, its first derivative with respect to 
|/3| 2 vanishes for \/3\ 2 = B* el, where 

B * _ 1 +P~ 2Vpcos(Ay?) 

4[p-^cos(A^)] ' [ ° j 

while it is strictly negative for \f3\ 2 < B* , and strictly 
positive for |/3| 2 > B*. Thus the fidelity F^ exhibits its 
minimum at \/3\ 2 = B* if B* G [0, 1] and at \(3\ 2 = 1 oth- 
erwise. The corresponding map for various combinations 
of p and Aip is shown in Fig. [TJa). Clearly, for a broad 
range of values for p and Aip, the fidelity is minimized for 
|/3 1 2 = 1, but it does also exist a non- negligible range of 
values for which the fidelity is minimized for input qubit 
states with |/3| 2 < 1. In the former case, the correspond- 
ing minimum value is F m i n (p,Aip) = p whereas in the 
latter 

- i + f COS<A rt - ^J^^y 

(17) 

Let us now compare F m - ln (p, Aip) to F^,(l,p, Atp) = p, 
which corresponds to the transfer of probability (excita- 
tion), as well as to F(p, Aip). The differences are shown 
in Fig. [Tfb) and Fig. Efc), respectively [28j]. Clearly, 
when the disorder is such that the random phase Aip is 
highly peaked around 0, we have F m i n (j>, Atp) = p, for all 
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FIG. 8. (Color online) Density plot of probability (|18l) as a 
function of the diagonal (a v ) and off-diagonal (cr^) disorder, 
for a chain of length (a) N = 12 and (b) JV = 21. The proba- 
bilities have been estimated for e = 10 -2 on 1000 realizations. 



the probability [H| 

Prob[p > 1/2 n |p - Av?)| < e] (18) 

for e <C 1, as a function of cr r; and <7£ for two dif- 
ferent values of N. As shown in Fig. [FJ in both 
cases the probability ([T8"|) exceeds 90% for weak disor- 
der, which means that in this case the transfer of the 
excitation (probability) can describe reliably the transfer 
of the state up to error e. By contrast, the probability 
Prob[p > 1/2 n \F - F mhl \ < 1(T 2 ], is practically neg- 
ligible and is not shown here. Once more therefore we 
see that F(p, Atp) is not a reliable measure for the QST 
under realistic conditions, since it tends to overestimate 
the performance of the protocol, relative to -Fining, A<p). 
V. CONCLUSIONS 



values of p e [0.5,1]. Hence, in this case one can safely 
consider the transfer of probability as a measure for the 
transfer of the state. For any other combination of Atp 
and p, we have F m - ln (p, Atp) < p, and thus F m i n (p, Atp) 
has to be used for the quantification of the QST in a sin- 
gle realization. By contrast, according to Fig. [TJc), there 
is only a very narrow regime of combinations for p and 
Ap, where \F min (p, Atp) - F(p, Atp>)\ < 0.02. This regime 
is actually very close to the ideal scenario, and is hard to 
be achieved in practice. 

It is worth emphasizing here that the main observa- 
tions of this section so far apply to any QST protocol 
for which the single-realization fidelity is given by Eq. 
(|T0)) . irrespective of the type of noise. In practice, for a 
certain QST Hamiltonian the parameters p and Atp vary 
from realization to realization, with p being always an 
upper bound on F min . Hence, to be on the safe side, 
one has always to quantify the performance of a QST 
Hamiltonian in terms of F m [ n . The above results, how- 
ever, suggest that it does exist a regime of diagonal and 
off-diagonal disorders, where F m i n seems to be well ap- 
proximated by p. Of course, the details of this regime 
depend on the QST Hamiltonian under consideration as 
well as on the modelling of the disorder. For instance, 
we have estimated numerically for the model of Sec. IIIII 



Wc have presented a statistical analysis on the perfor- 
mance of QST Hamiltonians in a single realization and 
in the presence of static disorder. It has been shown that 
the average-state fidelity, usually employed in the liter- 
ature for related studies, may fail to describe accurately 
the performance of QST Hamiltonians in a single realiza- 
tion. Most importantly, for a large class of input states 
it may also lead to faulty conclusions about the success 
of the transfer. For the sake of concreteness our analy- 
sis has been restricted to a particular well-studied QST 
Hamiltonian, and to a certain model of disorder. This al- 
lowed us to quantify the deviations of the average-state 
fidelity from the actual performance of the protocol in 
single realizations, as well as to identify the class of in- 
put qubit states and the size of the chains for which these 
deviations become fatal. 

Analogous conclusions are expected to be valid for 
other Hamiltonians and other models of disorder as well 
- albeit perhaps with some quantitative differences — 
and suggest that any studies on the robustness of QST 
Hamiltonians in the presence of noise have to rely on the 
minimum fidelity, while for weak disorder the transfer of 
the excitation (or probability) may be a rather reliable 
measure as well. In view of the present results, some 
of the previous studies on the classification of various 
QST Hamiltonians based on the average-state fidelity, 
may have to be revisited. 
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